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Abstract. By analyzing the dynamical properties of particle production in relativis- 
ts heavy ion collisions, it is possible to characterize the final stage of the equilibration 
process occurring in the collision fireball. In this work, we use the Hagedorn model 
coupled with non-extensive statistics to evaluate the transverse momentum spectra of 
positive pions for various event centrality classes in Au+Au collisions at RHIC with 
center-of-mass energies of 62.4GeV/A and 200GeV/A. We find that, by assuming an 
energy distribution that incorporates particle correlations, it is possible to explain 
the entire ir + pt spectrum as measured by RHIC. We find that spectra from central 
collisions, when compared to peripheral collisions, are consistent with a system that 
has smaller values of the non-extensivity parameter "q" and higher values of temper- 
ature. Comparison between different beam energies also shows a variation of the "q" 
parameter. The result is discussed using the interpretation that the "q" parameter 
is a measure of particle correlations within the system. Under this assumption, our 
results show that more central collisions are consistent with a system with less particle 
correlations. 

PACS Numbers: 25.75.-q; 25.75.-Ag; 24.10.Pa 



Introduction 



Relativistic Heavy Ion collisions have been studied for their potential in exploring 
many aspects of the QCD phase diagram. The particle production, as measured by 
experiments such as the Relativistic Heavy Ion Collider (RHIC) at the Brookhaven 
National Laboratory, exhibits transverse momentum spectra that are divided into 
three regions described by different physical models P]. For instance, particles with 
relatively low momenta (< 1.5GeV/c), which exhibit an exponential distribution, can 
be understood to be thermalized according to Boltzmann statistics. The intermediate 
Pt range, 1.5GeV < px < QGeV/c, might be considered as coming from quark 
coalescence and recombination, while the highest momenta pt > QGeV/ c could originate 
from parton fragmentation. Thus, transverse momentum spectra are usually explained 
differently according to the momentum range considered. 
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It would be of interest to find a single physical description for measured momentum 
spectra over the entire px range. This has been attempted before with different 
approaches. For instance, a framework of percolation of strings can be used to find 
a universality in the transverse momentum spectra [2]. Alternatively, temperature 
fluctuations have been shown to generate spectra compatible to the ones measured 
experimentally [3j. In this work, we chose to use a non-extensive statistical approach 
[HOE]. In this formalism the entropy is no longer additive. When computed in the 
canonical ensemble, the occupation probabilities are given by 

Pi = ^-[l-P (?-l)ei]5^ (1) 

with Z q being 

w 

Zq = J2[l-(3(q-l)e t ]— (2) 

i-1 

where Pi is the occupation probability of a state with energy e,, W is the number of 
accessible states, (3 is a real number and "g" is the non-extensivity parameter, so called 
because the composition rule for entropy deviates from additiveness by a mixed entropy 
term proportional to (1 — q). These expressions reduce to the Boltzmann occupation 
probabilities in the limit q — > 1, where entropy also becomes additive as in classical 
statistical theory. 

The probabilities given by equation [1] therefore give us an energy distribution, but 
it is still necessary to relate energy to transverse momentum. To do this, we assume 
a thermal equilibrium of a relativistic gas, as in the Hagedorn model, where energy is 
distributed according to equation[TJ One can then determine the momentum distribution 
to be p] 
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Here, T denotes a real number equal to 1//3 which is related to the slope of the 
resulting curve, pi is the longitudinal momentum and Pt is the transverse momentum 
of the produced particles, and /i 2 = Wq +Pt, where m is the rest mass, is the so called 
transverse mass of the produced particles. 

This function can be used numerically to adjust transverse momenta. It is, 
however, of crucial importance to understand the fit parameters in order to gain physical 
knowledge of the dynamics at hand. 



Dynamical Model 

A dynamical interpretation for the g-parameter has been suggested in reference [8] and 
proposes a discrete dynamical model which can be simulated easily. Consider a system 
consisting of ./V particles, each having energy E\ in multiples of some energy quantum, 
AE. The system evolves in discrete time steps in which two particles i and j are 
randomly chosen to be updated and one of the following two updates takes place: 
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Figure 1. a) Tsallis distribution fit to spectra generated using different r values in 
discrete dynamical model, b) Obtained q and T fit parameters according to the rvalue 
in the dynamical model. 



• With probability 1 — r, the two particles exchange one energy quantum, i.e. 
Ei = Ei ± AE and Ej = Ej =F AE. If the choice is such that a negative energy 
should arise, nothing happens and this step is skipped. 

• With probability r, the two particles are set to the energy of the lowest original 
energy, i.e. Ei = Ej = Min(Ei, Ej), and the excess energy is passed on to a third 
randomly chosen particle, effectively introducing a two-particle correlation. 

It is therefore possible, by changing the r parameter in the model, to increase or 
decrease the degree of correlation in the dynamics. We simulated such systems using 
N = 100 in a system saturated with energy, in which the average energy per particle 
is (E) = 180. In this case, at least 10 6 time steps are required in order to attain an 
energy distribution curve that is stationary It is noteworthy that different iV values 
change this relaxation time, as more iterations will be needed to attain equilibrium if 
more particles are present. 

This model produces energy distributions that are similar to Tsallis' energy 
distribution curves and can be fitted with equation HJ shown in Fig. UK The fit 
parameters indicate that the r value and the adjusted q value are closely related (Fig. 
[Tb) and the q value is indicative of a higher degree of correlation in the dynamics. Fig. 
[Tb also shows that the T and q parameters are anti-correlated, which can be seen as 
a reflection of conservation of energy. Populating higher energies with a small amount 
of particles implies taking this energy from the bulk of the other particles, effectively 
increasing the slope parameter and decreasing the T parameter. Other configurations 
varying iV and (E) were also studied and in all cases we observed the same qualitative 
behaviour for different r values. 

This toy model also makes the study of event-by-event energy fluctuations possible. 
By simulating many independent systems and looking at the energy occupation of the 
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Figure 2. Systematic study of event-by-event fluctuations, as measured by width 
of gaussian fits to energy occupancy using discrete dynamical model, a) shows data 
normalized by reference width and b) shows data normalized according to gaussian 
mean. The dashed line corresponds to the reference value for r = 0. 



2 — 15 energy quanta bin and fitting these data with a gaussian curve, it can be seen that 
the fluctuations increase by roughly 20 — 25% (Fig. [2^,) when increasing the correlation 
parameter r in the model. If we normalize the gaussian widths by the mean value, it 
can be seen (Fig. [2b) that this increase in the a values is not due to the central limit 
theorem for r > 10~ 2 , as this would require that the increase be proportional to the 
mean value. 

Comparing with Experimental Data 

We have looked at the 7r + momentum distribution as measured and published by the 
STAR experiment from pions from Au+Au and p+p collisions (see [9] and [10]) and 
fitted the Hagedorn gas model coupled with Tsallis' distribution, as in eq. El This can 
be readily compared to other full-range candidate functions, as seen in Fig. [3^. The 
best x 2 value obtained is from the Hagedorn+Tsallis model. 

This approach can be used systematically accross different centrality bins, as well 
as with p+p data, as shown in Fig. [3b. It is remarkable that the \ 2 values are good 
throughout all these datasets. Using data measured with a beam energy of 62.4GeV 
also results in the Hagedorn+Tsallis function having the best \ 2 values throughout all 
centrality bins studied. 

The dependence of the q and T parameters from these fits show that the more 
central collisions always exhibit a lower g-value when fitting, as well as a higher T 
parameter (Fig. HJ). The p+p data fits well in the systematic trend. Considering the 
Dynamical Model introduced in the previous section, where it was proposed that q could 
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Figure 3. a) Comparison of full-range fits to RHIC tt + data at 200GeV. b) Hagedorn 
model with Tsallis energy distribution fits accross several centrality bins and p+p data 
from RHIC at 200GeV. 
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Figure 4. q and T fit parameters obtained for 200GeV Au+Au, p+p and 62.4GeV 
Au+Au data, according to number of participants in the collision, calculated with 
optical glauber model. Larger number of participants correspond to more central 
collisions. 



be used as a measure of the correlation in the dynamics, the results shown in Fig. H] 
suggest that more central collisions have dynamics in which particles are less correlated. 
If the main source of correlation is assumed to be jet production, this could be an 
indication of jet quenching. The slight decrease in the T parameter is also consistent 
with an increased r value in the dynamical model. 

Comparing the fit parameters for the two different datasets measured with two 
beam energies, it can be seen that the higher beam energies exhibit a higher g-value as 
well as lower T parameters. Higher q values, using the same dynamical model analogy, 
translate to higher r value and a higher jet formation. 

It is noteworthy that this analysis has been done before ([?]) with e + e~ — > hadrons. 
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Figure 5. Comparison to fit results published for e + e — >■ hadrons by Bediaga et Al 
[7]. Only most central data from heavy ion collisions have been plotted. 



Our results can be compared to the previous as in Fig. El and even though the two 
physical processes are very different, in both cases the Hagedorn Gas Model coupled 
with Tsallis statistics is quite effective in describing the spectra obtained. The general 
trend of higher q- values and lower T-values is the same in both analyses. 



Conclusion 



We have shown that, by using the Hagedorn gas model with energies distributed 
according to Tsallis statistics, the full-range transverse momentum spectra from heavy 
ion and p+p collisions can be described. It is remarkable that this approach works well 
for both heavy ions and p+p, and throughout the entire measured range of transverse 
momenta. 

A preliminary interpretation of the fit parameters indicates that there are less 
correlations in the dynamics of more central collisions, as well as a higher slope parameter 
T. If we take our dynamical model and map its parameters to q and T, we see that this is 
consistent with a smaller r value. If we interpret correlations in the dynamics as coming 
mainly from jets, this indicates that in more central collisions there are proportionally 
less jets, i.e. jet quenching. A higher energy density has also been measured before for 
more central AuAu collisions at RHIC In addition, when comparing transverse 

momentum data from collisions at different beam energies, we see higher q parameters, 
which, when mapped to higher r values and using the same analogy, is an indication of 
dynamics in which more jets are produced. 

Further studies are warranted to fully understand the fit parameters as well as the 
physics that lies in the full-range transverse momentum spectra considered as such. It 
would also be of interest to investigate the spectra of other particle species using this 
same framework. Other frameworks have already been used to do this [12J, where the 



7 



blast wave model is used to describe various particle spectra. 

Finally, it is noteworthy that the transverse momentum range considered in this 
work is relatively large compared to the available transverse momentum data for other 
particle species due to experimental limitations. Our investigations show that having 
data all the way up to roughly 4 — 5GeV/ c is best for the parameters to converge ideally. 
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